Semiclassical spectral correlator in quasi 
one-dimensional systems 



Petr Braun 1 ' 2 , Sebastian Miiller 3 and Fritz Haake 1 

1 Fachbereich Physik, Universitat Duisburg-Essen, 47048 Duisburg, Germany 

2 Institute of Physics, Saint-Petersburg University, 198504 Saint-Petersburg, Russia 

3 Department of Mathematics, University of Bristol, Bristol BS8 1TW, United 
Kingdom 

E-mail: Petr.Braun@uni-due.de 
Abstract. 

We investigate the spectral statistics of chaotic quasi one dimensional systems such 
as long wires. To do so we represent the spectral correlation function R(e) through 
derivatives of a generating function and semiclassically approximate the latter in terms 
of periodic orbits. In contrast to previous work we obtain both non-oscillatory and 
oscillatory contributions to the correlation function. Both types of contributions are 
evaluated to leading order in 1/e for systems with and without time-reversal invariance. 
Our results agree with expressions from the theory of disordered systems. 

PACS numbers: 05.45.Mt, 73.21.-b, 73.20.Fz 
1. Introduction 

In the field of quantum chaos and disorder, the behavior of quasi one-dimensional 
systems such as long wires is clearly distinguished from the behavior of "normal" 
chaotic or disordered systems. Most importantly, quasi one-dimensional systems display 
Anderson localization PQ, i.e., wave functions are localized in only part of the wire and 
the conductance is suppressed. Anderson localization has important consequences for 
the statistics of energy levels [21 [3]. For normal systems the energy levels tend to 
repel each other; the spectral statistics is universal and agrees with predictions made by 
averaging over random-matrix caricatures of the possible Hamiltonians, according to the 
so called BGS conjecture [I]. In contrast the spectral statistics of quasi one- dimensional 
systems depends on the length (and thus the diffusion time To); in the limit of large 
length the energy levels belonging to the localized wave functions become independent 
and hence show Poissonian statistics. 

This difference between normal and quasi one-dimensional systems is well 
understood for disordered systems. Notable approaches are based on the DMPK 
equation [5J and on the nonlinear sigma model, a field-theoretical technique to evaluate 
averages over different realizations of the disorder potential. From the latter, localization 
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could be extracted in The appropriate definition of quasi one- dimensional behavior 

arising in this context is that the classical diffusion time becomes comparable to or 
larger than the relevant quantum time scales, in particular the Heisenberg time Th = ^ 
where A is the mean level spacing. A random-matrix model for systems of this type 
was considered in [HI El HQ] . 

For clean chaotic systems (e.g. wires in which the classical motion becomes chaotic 
due to the shape of the boundary) the effects of quasi one-dimensionality are less well 
understood, and most of the literature is restricted to normal systems. A quantity that 
has attracted a lot of attention in this context is the spectral correlation function R(e). 
For dynamical systems with a level density p(E) this correlation function is defined by 

*< e ) = A><„(* + (!) 

where e is a real energy offset, the brackets denote an average over the center energy 
E and A is the mean level spacing. For simplicity we assume that A is brought to 1 
by appropriately scaling the energy levels. Random matrix theory (RMT) now makes 
predictions for R(e): For systems without time reversal invariance an average over the 
Gaussian Unitary Ensemble (GUE) of RMT gives — 7^2 + while for time reversal 
invariant systems an average over the Gaussian Orthogonal Ensemble (GOE) leads to 
infinite power series in \ and cos „ 2e . The slow (power law) decay of oscillations leads to 
a singularity in the Fourier transform of R(e) (the spectral form factor) at time t = Th- 

To show that individual systems are faithful to these predictions a semiclassical 
approach was proposed in [UJ, [12j [T3] . The essential idea is to express p(E) as a sum 
over periodic orbits, using Gutzwiller's formula [T4], and then study the interference 
between contributions of these orbits. The leading non-oscillatory contribution arises 
from "diagonal" pairs of identical (up to time reversal) orbits [UJ. The remaining terms 
were accessed only recently in [13], [151 EES EH] ■ 

In the present paper we want to generalize these new results to quasi one- 
dimensional systems. For these systems the diagonal approximation to the small-time 
form factor (the Fourier image of the non-oscillatory part of R(e)) was evaluated by 
Dittrich [TH]. First results on off-diagonal contributions were obtained by Schanz and 
Smilansky [19] for one-dimensional quantum graphs, and by Brouwer and Altland [20] 
who semiclassically explained localization for quasi one-dimensional systems modelled 
by an array of quantum dots. In contrast, we will focus on general quasi one- dimensional 
systems such as long wires. We use a periodic-orbit expansion not of the correlation 
function itself, but of a generating function which yields R(e) upon taking derivatives 
[T6| [TT] . This enables us to determine, to leading order in -, both the non-oscillatory 
and the oscillatory parts of R(e). In this order we see that the effects of quasi one- 
dimensionality reduce to modification of the periodic orbit sum rule suggested in 
[18J. For systems without time reversal invariance it suffices to perform a diagonal 
approximation on the level of the generating function. In contrast, for time- reversal 
invariant systems this diagonal approximation still captures only the non-oscillatory 
part; the evaluation of the oscillatory part involves off-diagonal contributions of pairs 
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of non-identical but similar orbits. Both for systems with and without time-reversal 
invariance we reach agreement with results for disordered systems by Andreev and 
Altshuler [3]. Our results illustrate how semiclassical methods are useful not only for 
describing universal features of "normal" systems but also deviations from universality. 

Higher-order corrections in 1/e should be similarly accessible; their calculation 
needs taking into account more complicated groups of correlated orbits introduced in 
previous work on normal systems [T3| 115], combined with treatment of higher order effects 
of quasi one- dimensionality. An extension of this approach to Anderson localization 
appears within reach. 



2. Two-point spectral correlator and the generating function 

To get started we briefly review how the correlation function R(e) can be accessed 
through a generating function. Following Ref. [33] we write R(e) as the real part of the 
complex correlation function 

e ± =e±i 7 , (2) 
R(e) = limReC(e+) 

and determine the latter from a generating function, the energy-averaged combination 
of four spectral determinants 

V A B C D) \det (E + e+- H) det (E + e B - H) / K J 



as 

, . 1 Pi* 7. 

(4) 



1 a Z 

R(e) = lim Re C (e + ) = — + 2 Re lim 



2 7-0 de\de B 

Here the subscripts ± indicate small positive or negative imaginary parts. The symbols 
||, x denote two alternative ways of identifying the energy arguments, to be referred to 
as "columnwise" (||) and "crosswise" (x), 

II : € a = e c = e A' e B = e D = ~ e+ columnwise , (5) 

x : e\ = e + , e B = — e + , = — e~, = e _ , 7 — > +0 crosswise . (6) 

Both procedures would yield the same result for the two-point correlator if 
implemented rigorously. However, we shall have to calculate Z semiclassically, and that 
approximation entails two different expressions, one (||) reproducing the non-oscillatory 
part and the other (x) the oscillatory part of R(e). To obtain the full result both 
expressions have to be added. In jTTj it was shown that this addition can be understood 
naturally in terms of an improved semiclassical approximation preserving the unitarity 
of the time evolution (the Riemann-Siegel lookalike formula [2T]). 
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The semiclassical approximation for Z is based on Gutzwiller's formula for the trace 
of the resolvent 

tr(£+ - H)- 1 = + J2F a e lS ^ h (7) 

a 

The factor proportional to E + is the smooth (Weyl) part of the level density and the sum 
taken over periodic orbits with S a , T a , F a action, period and stability coefficient of the 
ath orbit. Integration then yields the semiclassical approximation of the determinant 

det (£ + - = exp [ - j E alETi (E - H)' 1 ' 



exp 



v A 

Substituting such expansions for all four determinants in Z and expanding, e.g., 
S a (E + e%) « S a (E) + T a e\ we obtain 

T a + ■ T a C + N 



e m-*B-4+*n) exp [ £ Fa e k&W - e^ e ° ) (8) 



3. Systems without time reversal invariance 

3.1. Diagonal approximation 

The semiclassical representation (jSj) falls into a product over periodic orbits, 

Z = e K4-^-4+^) Zo j Zq = J[z a , (9) 

a 

z a = exp [F a (e^ t+A - e^ e ^ e ^ Sa{E) 



- f 

-J A 



AC 



+ F* [e~' l ^ e ~ B -e~'^ eD ) e -^ Sa{E) ] 



— fa* 
J BD 

The diagonal approximation [11] assumes that for systems without time-reversal 
invariance contributions of different periodic orbits are uncorrelated. For the generating 
function this means that the energy-averaged Z becomes a product of single-orbit 

averages, 

(^o)diag = II<*> • ( 10 ) 

a 

The energy average in (z a ) does away with rapid oscillations due to the phase 
factors e ±lS ( E " h in the exponent. Since z a is periodic in the phase $ a = S a (E)/h we 
may average with respect to $ a , over a single period 2tt. This yields 

{za) = h r ^ exp ( ; ^ ei * + / ^ e_i$ ) = j ° (v^^d) , in) 
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where Iq is the imaginary-argument Bessel function. The expansion In I [y) = ^ 

64 i •••b"™ \~u/ — I L^a J AUJ BU 4 



o J + . . . gives (Z ) = exp [£ a f a Ac f%* D - \ EaifXcfBb) 2 + •••]• In the semiclassical limit 
(Tff — *■ 00) it suffices to keep only the leading quadratic term in the exponent||| 
Our task is thus reduced to calculating the periodic-orbit sum 

ln(^o) diag = E/Ac/BD (12) 

= El^l 



I 2 ( e 1 T H € A _ e l T H £ C 



3.2. Diffusive vs ergodic behavior 

For "normal" systems the sum over periodic orbits can be done using the well-known 
sum rule of Hannay and Ozorio de Almeida [22] 

poo dT 

El^l 2 (-) = / i-(-) (13) 

which expresses the approximately ergodic behavior of long orbits; short orbits below a 
certain classical period T have to be excluded. 

The quasi one-dimensional character of long wires requires a modification of 
that sum rule [18J. The momentum components in all directions and the transverse 
coordinates can still effectively randomize after a few bounces against the boundary, 
but the longitudinal coordinate takes much longer to explore the whole length L of the 
wire. Given chaos inducing boundaries, each bounce will with finite probability change 
the sign of the longitudinal momentum component and thus entail diffusion along the 
wire. The longitudinal coordinate becomes randomized only for times long enough to 
explore the wire, i.e., times in excess of the diffusion (Thouless) time To = L 2 /D 
where D denotes the diffusion constant. Orbits with periods T much smaller than that 
Thouless time will have explored only the small fraction y/DT/L of the whole length, 
and the inverse of that fraction must be expected as a factor of increase of the right 
hand side of the above sum rule, relative to orbits with periods T 3> Tjj. 

Dittrich [18] has determined the aforementioned factor of increase for arbitrary 
values of the ratio T/T D as the integral J L dx p xo (x , T\x ) = P(T) of the probability 
density of return to an arbitrary point Xq of departure after a time T, for a one- 
dimensional random walk. Solution of the diffusion equation for a wire of length L, 
led to the "enhanced-return" factor 

P(T) = E exp ( - = \[l + ^3 (0, e"^)] (14) 

n=0 Z1 D 7 i L 

| For this term the decrease of the stability coefficients \F a | 2 with the period is just compensated by the 
exponential increase in the number of orbits. For all other terms F a decreases faster; the contributions 
of the long orbits then become exponentially small whereas the shortest orbits make for corrections of 
the order o(T /T H ). 
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where ^(u, 9) denotes the elliptic theta-function of the third kind [23J. The sum rule 
modified for quasi one dimensional systems thus reads 

ei^i 2 (-)= r^p{T){-). as) 

a JT 1 

Hannay's and Ozorio de Almeida's sum rule is restored if To is so small compared 
with T that only the n = term in ffl4|) survives and P(T) ~ 1. In the opposite 
limit T D > T we have P{T) ~ ^T d /2ttT = L/y/2-irDT which agrees with the above 
qualitative expectation. 

Inasmuch as orbit periods of the order of the Heisenberg time T# determine the 
level statistics the time scale ratio 

^ < 16 > 

will play an important role for the wires under study here. In particular, the crossover 
from normal behavior (£ 3> 1) to the quasi one-dimensional behavior under study here 
takes place for ( of the order unity; for ( <C 1 the spectral statistics must be compatible 
with localization. 

We now invoke the modified sum rule and the identity / °° ^(e mT — e lbT ) = ln|; 
the lower limit T of the time integral could be replaced by zero, accepting a negligible 
error o(T /T H ). We thus get the diagonal approximation for our generating function as 



^ f°° dl 

n=CT U 1 



°° dT -nlJfT 



+ e _i 2^-( £ D-4) _ e -i^( £ D-<4) 



= X>f^ ~ — ^vr 1 ~ — rf ( 1T ) 

„=0 (lC^ 2 + 4 - e B) [< n2 + e C ~ e D) 

The infinite product implicitly involved here can be brought to a closed form using 

ft ~2~7T = ^717' = v^sinhTrv/a 7 , 

^ n 2 + b if [b) 

whereupon we arrive at our final result for the generating function in the diagonal 
approximation, 



(Z) diag = eK4-^-^B) 



iC J r y iC 

In the limit Tp — > or £ — > 00 we have <p(x) — * 7rx, and the generating function tends 
to its familiar form for normal systems [T5] . 

5. 5. Two-point correlator and form factor 

Substituting (Z) diSLg for (Z) in (TjJ and identifying energies columnwise (||) we obtain 

00 1 

C11 (e) = - 2 V 9 

11 (i(n 2 + 2ef 
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1/1 1„ „ ieTT 2 \ , s 

= o - + -0 cot 9 + — 5- , 19 

2e 2 \2 4 2Csin 2 ^y ' v ; 

9 = {l + i)*Jl. 

Upon taking the real part we are led to a cumbersome expression for R(e) equivalent 
to the earlier RMT [3] and semiclassical [18] results for the non-oscillatory part of the 
correlator. In the limit ( — > oo the GUE behavior i? non _ osc (e) = — l/2e 2 is restored. 
The crosswise (x) identification of parameters, on the other hand, entails 

2vr 2 e i2e 

C y (e) = — r — — , (20) 



COsh 27TW^ — COS 2lTy 



and now the real part yields the oscillatory part of the correlator R osc (e), in agreement 
with what Andreev and Altshuler {3] had found through an average over an ensemble 
of disordered systems. The GUE expression R osc (e) = cos2e/2e 2 follows in the limit 
C -> oo. 

For finite ( the amplitude of oscillations of the spectral correlation function tends 
to zero exponentially with e — > oo instead of the power law characteristic of normal 
systems. That leads to qualitative changes in the spectral form factor K{r) where r is 
the dimensionless time, r = t/Tjj. For r > 0, the form factor can be defined as the 
Fourier transform, 

K{t) = — e- i2£T C{e)de. (21) 



27T J — 00+17 

In normal systems without time reversal invariance K(r) experiences discontinuity of 
its first derivative at r = 1 introduced by the Fourier transform of the oscillatory part 
of the spectral correlation function. In quasi one-dimensional system this discontinuity 
is replaced by a smooth transition from the small-time to large-time behavior. Another 
change associated with the non-oscillatory part of the spectral correlation function is the 
much faster growth of K(t) at small r: The respective ("parallel") part of the diagonal 
form factor first deduced semiclassically in 



K ll (r)=rJ2e- (nr = rP(T H r) (22) 

n=0 

grows like square root rather than linearly. That faster rise toward the saturation value 
unity may be seen as a (slightly indirect) hint to localization; in the Poissonian limit 
this rise would be a jump, K(t) = 1 for all r > 0. 

As a note of caution it has to be mentioned that, in contrast to normal systems 
without time reversal, the diagonal approximation for the correlation functions no longer 
coincides with the exact result in quasi one-dimensional systems. In particular, the 
Fourier transform of the oscillatory part (120]) is no longer zero for r < 1 tending to a 
finite negative value for r — ► +0. Consequently the total form factor of the diagonal 
approximation becomes negative for small r although the exact form factor is known to 
be non-negative. 
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4. Time reversal invariance 

Interesting changes arise if time-reversal invariance holds. Periodic orbits then exist in 
time reversed pairs (a, a) with exactly the same action and stability coefficients. The 
generating function becomes a product of contributions of different pairs, and these 
pairs are uncorrelated in the diagonal approximation. Due to z a = z a the generating 
function is squared compared to ({TBI , the Weyl factor apart, 



<P V 



Using this generating function with columnwise identification of arguments we find that 
the non-oscillatory part of the two-point correlation function and the small-time form 
factor are doubled compared to ( fT9l) and (122]) : this is in line with Refs. [3] and [TB] . 

Remarkably, for time-reversal invariant systems the diagonal approximation yields 
no oscillatory contributions to the correlation function, i.e., there are no terms of order 
cos 2 2e . This can be understood as follows. In the crosswise limit ([6]) we have 

et ; -e- B = e + A -e- D = \2 1 , 7 ^0 (24) 

such that we can replace <p(x) — > nx in the numerator of (Z) dia _ g ; this gives 

(Z> diag oc(e+-e B ) 2 ( e +-e D ) 2 (25) 

such that (Z)diag tends to zero like 0(7 4 ). The two derivatives w.r.t. e^,e B can only 
eliminate two factors 7. This leaves a result that tends to zero like 0( 7 2 ) which means 
that Cdiag,x = 0. 

To derive the oscillatory component of the spectral correlator we thus have to 
go beyond the diagonal approximation and take into account correlations between the 
factors z a in the generating function related to different periodic orbits. For the relevant 
correlated orbits the differences (z a z b ) — (z a ) (z b ) , (z a z b z c ) — (z a ) (z b ) (z c ) etc. must 
be non-zero. In view of the semiclassical limit this is possible only if the respective 
actions have a chance to cancel, e.g. if S a (E) « S b (E) or S a ~ S b (E) + S C (E). 
Such correlations between orbits indeed exist for chaotic dynamics; they stem from 
"encounters", i.e., places where two or more stretches of the same orbit or of different 
orbits are close and almost parallel or antiparallel to each other. By changing the 
connections inside these encounters one can turn, e.g., an orbit a into an orbit b with 
almost the same action, or split it into two orbits b and c whose sum of actions is 
close to the action of a. We shall refer to such sets of correlated orbits as "bunches". 
The simplest encounter involves two almost antiparallel orbit stretches; the bunch it 
generates is the famous Sieber-Richter pair (containing one orbit where the encounter 
forms a crossing in configuration space and one where it forms an avoided crossing) [13J. 
More complicated scenarios were introduced in [151 [IB"] . It has been shown in [16] that 
taking into account both the "diagonal" correlations and those related to bunches gives 
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correct semiclassical asymptotics of the generating function and the correlation function 
of normal systems. The generating function was found as 

(Z) = (^)diag (1 + (ZU) , (26) 

where the off-diagonal part (Z) Q g contains the contributions of the bunches mentioned. 

Let us now determine the term in (Z) g responsible for the leading oscillatory 
contribution to the correlator. Multiplication with this term must remove one factor 
£c~ € b ~ > ^7 and one factor e\ — — ► 227 from (1251) . The product is then proportional 
to (e\ — e^)(ej — e#) and survives differentiation w.r.t. e\ and and taking the limit 
7 — > 0. For normal systems the term of lowest order in - satisfying this condition reads 
(see Eqs. (13,14) of the on-line version of Ref. [16J) 

(Z) oS , x = ~ (+ * y (27) 

In a quasi one-dimensional system, due to the diffusive dynamics, Eq. ( 12T1) would 
be replaced by an expression analogous to the diagonal approximation, i.e. a summand 
as in (T2"7j) . plus terms where the energy differences in the denominator are shifted by 
finite imaginary amounts of the type i(n 2 , with integer n as in ( 1T7I) . However, the shifted 
terms would, for n 7^ 0, no longer diverge in the limit (x). Hence combined with the 
above factor 0(7 2 ) they would yield vanishing contributions to the correlation function. 
Therefore the leading term in C x (e) is still due to (1271) . Substituting (1271) into (1261) and 
calculating derivatives in the crosswise procedure, we obtain the oscillatory component 
of the complex correlator, 

87r 4 e i2£ 

lim C x (e) = ^— ? . (28) 

e 2 ( 2 ( cosh 27r,M — cos 2tt J 

Its real part coincides with the RMT result for the two-point spectral correlation 
function in the presence of time reversal[3], now deduced semiclassically for individual 
chaotic quasi one- dimensional systems. In the normal-system limit ( — > 00 (f28l) tends 
to the random-matrix expression exp (i2e) /2e 4 . On the other hand, for all finite ( 
the amplitude of oscillations diminishes exponentially with the growth of e. As a 
consequence the discontinuity at r = 1 of the third derivative of the GOE spectral 
form factor K(t) [9] is smoothed out in the quasi one- dimensional case. 
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